The flight of the bumblebee: 
vacuum solutions of a gravity model with vector-induced 
spontaneous Lorentz symmetry breaking 
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We study the vacuum solutions of a gravity model where Lorentz symmetry is spontaneously 
broken once a vector field acquires a vacuum expectation value. Results are presented for the 
purely radial Lorentz symmetry breaking (LSB), radial/temporal LSB and axial/temporal LSB. 
The purely radial LSB result corresponds to new black hole solutions. When possible, Parametrized 
Post-Newtonian (PPN) parameters are computed and observational boundaries used to constrain 
the Lorentz symmetry breaking scale. 
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I. INTRODUCTION 

Lorentz invariance is possibly one of the most funda- 
mental symmetries of Nature. It is theoretically sound, 
it has been extensively tested, and so far no clear cut 
experimental evidence has emerged to bridge its validity 
(see e.g. Refs. [1, 2]). It is therefore no surprise that 
most theories of gravity encompass this symmetry, and 
little attention has been paid in understanding the im- 
plications of the breaking of Lorentz invariance in this 
context. 

There is, however, a window of opportunity for fiddling 
with a less stringent approach: as for gauge symmetries 
in field theory, many relevant after effects can arise if 
one assumes a spontaneous breaking of Lorentz invari- 
ance [3, 4]. In the context of gravity, the breaking of 
Lorentz invariance can be implemented if a vector field 
ruled by a potential exhibiting a minimum rolls to its vac- 
uum expectation value (vev), in the fashion of the Higgs 
mechanism [5]. This "bumblebee" vector thus acquires 
an explicit (four-dimensional) orientation, and preferred- 
frame effects may arise. 

The possibility of violation of the this fundamental 
symmetry has been widely discussed in the recent lit- 
erature (see Refs. [1, 2]). Indeed, string/Mtheory allow 
for a spontaneous breaking of Lorentz symmetry, due to 
the existence of non-trivial solutions in string field theory 
[3, 4], loop quantum gravity [6, 7], quantum gravity in- 
spired spacetime foam scenarios [8] and noncommutative 
field theories [9, 10]. Also, LSB could result from space- 
time variation of fundamental coupling constants [11]. 
Experimental tests of this symmetry breaking could be 
achieved, for instance, in ultra-high energy cosmic rays 
[12]. 

Efforts to quantify an hypothetical breaking of Lorentz 
invariance have been mainly concerned with the phe- 
nomenology of such spontaneous Lorentz symmetry 
breaking (LSB) in particle physics. Only recently its 
implications for gravity have been more thoroughly ex- 
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plored [5, 13]. In that work, the framework for the LSB 
gravity model is set up, developing the action and using 
the vielbein formalism. 

In this study, we focus on the consequences of such 
scenario. For this, one assumes three relevant cases: the 
bumblebee field acquiring a purely radial vev, a mixed 
radial and temporal vev and a mixed axial and temporal 
vev. For prompt comparison with experimental tests, we 
shall write our results in terms of the PPN parameters, 
when possible. 

The action of the bumblebee model is written as 



S = 




(R + £B»B v R fiU ) 
ViB^B^ ± b 2 )} 



(1) 



where k — 8irG, £ is a real coupling constant and b 2 is a 
real positive constant. The potential V driving Lorentz 
and/or CPT violation can be chosen to have a minimum 
0. 



at B^B^ ± V 



II. PURELY RADIAL LSB 

In this section we develop a method to obtain the exact 
solution for the purely radial LSB. For this, we assume 
a static, spherically symmetric spacetime, with a met- 
ric Qjiv — diag(—e 2 ^, e 2p , r 2 , r 2 sin 2 8), where <f> and p are 
functions of r. Also, we admit that the vector field B^ 
has relaxed onto its expectation value b^. It is imposed 
that D^by = 0. This enables us to calculate 6 M , using 
the affine connection derived from the metric g^ v ; for 
this, and since the only non trivial covariant derivative is 
with respect to the radial coordinate, it is assumed that 
6 M = (0, b(r), 0, 0). Hence, from 



P Q b 



(2) 



2 



follows that 



Hr) = C 1/2 b e p 



(3) 



where the is introduced for later convenience. It 
can be immediately understood that, as expected, b 2 = 
b^b^ = b 2 ,^ 1 is constant. The action can be thus written 
as 



R_ 
~2~k 



+ (g rr ) 2 b 2 (r)R r 



(4) 



3-L = 



3 + 6g^, b 2 



x x 2 



~ S — 

'-2 ~ 2 



(11) 



and hence L ~ 6§/2. From the definition of \f, we obtain 



e 2 "= (l-*or- 1+i ) _1 . 



(12) 



Comparing with the usual Schwarzschild metric, one 
can write 



where the subscript s stands for the spatial quantities. 
With the assumed metric, the determinant is given by 
y/—g = r 2 e p+ ^\ the scalar curvature and the relevant 
Ricci tensor non-vanishing component are given by 



2G L m L 

9rr = 1 r 



.9oo 



i+ 2G^n rL 



(13) 



R 



2 [1 + (2rp' - l)e~ 2 P] 



(.) 



where the prime stands for derivative with respect to 
r and we have integrated over the angular dependence. 
Also, 



where Gl has dimensions [Gl] = L (in natural units, 
where c = K = 1). One can define Gl = Gr^ L , where r 
is an arbitrary distance. The L — > limit yields G^ — > G 
and the usual geometrical mass, M = Gm, with dimen- 
sions of lenght. From now on we express all results in 
terms of M. The event horizon condition is given by 



((1')^ = ^ , (6) 

where b r is the (contravariant) radial component of 
bfj,. We now introduce the field redefinition = 

(1 - e~ 2 P) r- 2 , so that 



2M r L 

.9oo = -1 + =0 , (14) 

r r 

thus r s = (2Mr ~ L ) 1 / 1 ~ L . One can compute the scalar 
invariant I = R^u P \R pupX (the norm-square of the Rie- 
mann tensor), obtaining 



2p'e 



-2p 



2* 

r 



(7) 



and thus 



2p'e 



= 2* + r*' , 



R = 2(3* + rtf') , 
(& r ) 2 i? rr = 6 2 (2* + r*') . 



(8) 



We can now work with the action (see [16] and refer- 
ences therein): 



dre p+4 >r 2 



(3 + b 2 )^ + (l + ^)r^ 



(9) 



Variation with respect to <j) yields the equation of mo- 
tion 



(3 + 6 2 )* + (+l + ^)r*' = , (10) 

which admits the solution \I>(r) = * r _3+i , where we 
define 



48 



17 
12 J 



-L 6 



12 



M 



r 
Jo, 



2L 



* I1--L 



2L 



Io 



(15) 



where Io — 48 M r is the usual scalar invariant in the 
limit L — > 0. Hence, since J(r = r s ) is finite, one sees 
that the singularity at r — r s is removable. By the same 
token, the singularity at r — is intrinsic, as given by 
the divergence of the scalar invariant. Thus, an axial 
LSB gravity model admits new black hole solutions whose 
singularity is well protected within a horizon of radius r s . 
One expects a Hawking temperature given by 



T - 



1 



kf> Airr. 



= (2Mr L )- L/{1 - L) T ~ (2Mr " 



~ L Tn 



(16) 

where T = h/8nkBM is the usual Hawking temperature, 
which is recovered in the limit L — > 0. 

Of course, this exact solution does not allow for a PPN 
expansion, as the obtained metric cannot be expanded 
in powers of U = M/r. One can, however, compare 
with results for deviations from Newtonian gravity [17], 
usually stated in terms of a Yukawa potential of the form 



3 



V Y (r) 



Gym 



(17) 



Equating this potential to the one arising from <?oo in Eq. 
(13), one gets 

G L r L =G Y (l + ae- r / x ^j . (18) 
Expanding to first order around r — r yields 



G L r% (l + L^J =G Y (l + a-a^ , 



(19) 



which allows us to identify A = r and a = — L (with 
GV(1 - L) = G L r£ =G). This states that the effects of 
a radial LSB, probed at a distance ro from the source, 
can be interpreted as due to a Yukawa potential of cou- 
pling strength a = — L and range A = r . The negative 
sign of the coupling shows that the radial LSB yields a 
"repulsive" component (for r > ro), as can be seen from 
Eq. (13). 

Notice that this identification of the LSB effect with 
a Yukawa potential constraints the length A to be equal 
to the distance ro at which one tests deviations from the 
inverse square law. This is not the case with a "true" 
Yukawa perturbation, where each test of gravity at a dis- 
tance ro yields an allowed range for both a and A (al- 
though, of course, the test is sensible to deviations only 
at a scale A close to r . Hence, to obtain a bound on 
L one must look at the allowed value of a for the fixed 
A = ro at which an experiment is carried out. The most 
stringent bound is derived from planetary tests to Ke- 
pler's law, with Venus yielding A = r = 0.723 AU and 
L = \a\ < 2 x 1CT 9 [17]. 



III. RADIAL/TEMPORAL LSB 

We consider now the mixed radial and temporal 
Lorentz symmetry breaking. As before, we assume 
that the bumblebee field has relaxed to its vac- 
uum expectation value. Assuming temporal variations 
to be of order of the age of the Universe Hq , where 
Hq is the Hubble constant, one can, as before, con- 
sider a Birkhoff static, radially symmetric metric g^ — 
diag(—e 2 ^ 7 e 2p ,r 2 ,r 2 sin 2 (9)). Imposing as (physical) 
gauge choice the vanishing of the covariant derivative 
of the field B^, one gets b r (r) = £~ 1 / 2 A r e p and, sim- 
ilarly, &o(r) = £,^ 1 ^ 2 Aoe^, with Aq and A r dimension- 
less constants. One immediately sees that, as expected, 
b 2 = b p b^ = (A 2 - Al)^- 1 is constant. 

Since one now has both a radial and a time component 
for the vector field vev, the symmetry 6 — —p does not 
hold. Therefore, one cannot use the previous spatial ac- 
tion formalism depicted in Eq. (9). Instead, one resorts 
to the full Einstein equations, 



Gp. v — £ 



-{b a ) 2 R aa g — b^b v R^ v — b v b p R vll _ 



• (20) 



The additional equation of motion for the vector field 
vanishes, since it has relaxed to its vev and therefore 
both the field strength and the potential term are null. 
Introducing the metric Ansatz and the expressions for b^ 
one gets, after a little algebra, 



Goo 



3AlR 00 - A 2 r e 2 ^-P^R r 
A 2 e 2{p -^R 00 - 3A 2 r R r 



(21) 



We now write G o = 3o( r )e 2 ^ p \ G rr = g r {r), R 00 
/ (r)e 2 ^-' 5 ) and R rr = f r (r), where 



Mr) ee 




fr(r) EE 


(2 + rd')p> 


r 


9o (r) ee 


-l + e 2p 2p' 


~2 + V 


9r{r) EE 


1 - e 2p 26' 


~^2— + ^r 



(22) 



Inserting the above into the Einstein equations, it follows 
that 



9o(r) - -[3A 2 f (r)-A 2 r f r (r)] , 

9r(r) = \[Alf (r)-2,A 2 Mr)\ . (23) 

Hence, one must solve this set of coupled second order 
differential equations, with boundary conditions given by 
6{<x>) = p(oo) = 4>'(oo) = p'(oo) = 0. 

Before continuing, we point out that LSB is clearly 
exhibited: noticing that go + g r = fo + / r , one has 
(1 — 2^4o)./o = — (1 + 2A 2 .)f r , which is an explicit man- 
ifestation of the breaking of Lorentz symmetry. In the 
unperturbed case Aq = A r = 0, fo = —f r and one recov- 
ers the Schwarzschild solution 6 = —p from the equation 
.9o + 9r = 0. This symmetry does not hold in the per- 
turbed case, which yields fo ~ —(1 + 2A\ + 2A 2 )f r . 

To solve the set of coupled differential equations Eqs. 
(23) one considers an expansion of the metric in terms of 
6 — 6o + 56 and p = —6o — 5p, where 6o is given by the 
usual Szcharzschild metric: 



(24) 



and Sp, 56 are assumed to be small perturbations. Hence, 
one gets to first order 



4 



M 



fo(r) = -6<f>' 



^ _ 2M 



(3<5<// + 8p') + 5<t>" 



(25) 



and 



5p> 



1-4A 



A + B 



1 + 45 r 2(1 + 4B) 
which, after integration, yields 



(32) 



M 



/ r (r) = --V- I fWW + V)-^" • (26) 

r 

As expected, the above quantities are homogeneous in 
<5/9, S(f> and their derivatives. 

For the calculus of go(r) and g r (r), one first computes 
the contribution of the exponential term: 



1 _ e 2 P I _ e -2(4>o+S P ) 



1 - e - 2 ^(l - 2(5/?) 



(27) 



2Af 



1 - 



2M 



+ 



2 <5p 

^ 2M ^2 



Thus, one finds 



9o(r) 



and 



«5p 



0r(r) 



^ _ 2M r 



1 8p 



+ 6p' 



(28) 



r 1- 



2M 



(29) 



which are also homogeneous. 

To solve Eqs. (23), one first obtains a relation between 
5p and 8<j). For this, one sums both coupled differential 
equations, which yields 



-(54)' - Sp r ) = M-(A5<t>' + BSp') 
r r 



(30) 



M 



+ (A + B) 



1 - 



2M 



(3,50' + 8p') + (A + B)S<f>"} , 



where the simplifying notation A = A 2 ,, B = A 2 has 
been used (A and B are dimcnsionless). Since the LSB is 
presumably a small effect, one expects A <C 1 and B <C 1. 
Hence, one has 



5p> 



1 + 45- 



1-4A- 



A + B f 

2 1 - 2M 

r 

3(A + 5) f 



(31) 



1 - 



2M 



A + 5 



Dropping terms smaller than O(A), O(B), this expres- 
sion simplifies to 



2 + B-7A :± A + B 



2(1+45) 



2(1 + 45) 



(33) 



One can now introduce this expression in one of the 
Eqs. (23) and obtain, after a somewhat tedious compu- 
tation, the following ordinary differential equation 



-C x r 2 S$' + C 2 rScp' + C?,8<t) = 



(34) 



with 



Ci = A + 35 + AB + 95 2 ~ A + 35 
C 2 = 2 + 5-3A+16AB-2 , 
C 3 = 2 + 5-7A-2 . 

This equation has the solution 



(35) 



C 1 +C 2 -y/(C 1 +C 2 ) 2 +4C 1 C 3 

8(f) = K-r 2C i 



C 1 +C 2 + y / (C 1 +C 2 ) 2 +4.C 1 C 3 

+ K + r 2C i 



(36) 



where K_ and K + are integration constants. Since all Ci 
are positive, one obtains that y/ {C[ + C'2) 2 + AC1C3 > 
C\ + C2, and hence the second power-law term diverges 
when r — > 00. Since one demands that 8<p — ^ 0, this 
implies that K + = 0. The remaining first power-law 
term automatically satisfies the conditions <50 — > and 
8(f)' — > for all K_ = K. Hence, the perturbation is 
simply given by 8(f) = Kr~ a , where 



-C 1 -C 2 + y/(C 1 +C 2 ) 2 +AC 1 C; 
a = > . (37) 

and K has dimensions of [K] = L a . One can linearize 
this exponent with respect to C\ <C 1, that is 



C 3 



2 - 7 A + B 



d+C 2 2 + A + 3B ' 
so that a ~ 1. One can now compute (5p, obtaining 



(38) 



<5p = 



2 + 5- 7,4 A + 5 

2(1+45) ^ 2(1 + 45) 1 

l + «Cd±*>Ur- 



(39) 



Hence, the non-trivial components of the metric read 



■5 



9tt 



= 2(0 o +«50) 



-2(</> +<5p) 



,2Kr~ 



-{2+a{A+B))Kr~ 
1 2M 



-2K r 



I — 2M 



where one defines K r = [1 + a(A + B)/2]K ~[1 + (A + 
B)/2]K. To compute the PPN parameters, one first per- 
forms a Lorentz transformation to a isotropic coordinate 
system, that is, to one on which all spatial metric com- 
ponents are equal. Since the angular coordinates are not 
affected by the LSB dynamics, this amounts only to a 
change of the radial parameter, r — > £ = Thus, 
instead of explicitly deriving the associated metric g^ v 
through its transformation properties, one resorts to the 
invariance of the interval ds 2 , obtaining 



9tt(r)dt 2 +g rr {r)dr 2 + r W 

= gtt&dt 2 + xmde + edn 2 ) 



(41) 



Equating the angular part, one obtains A(£) = r 2 /£ 2 . 
Since the temporal part is equal on both sides of the 
equation (merely expressed in terms of r or £), one also 
obtains g rr dr 2 — A(£)d£ 2 and 



gli 2 d£, 

dr = — 

r £ 



(42) 



implying that 



/ 



-K r r~ 



dr = log £ + const . 



(43) 



dr = 



Co 



-rfCo • (46) 



Expanding the integrand around a = 1 (through r a — 
r + (a — l)r log(r)) and K r — 0, it can be easily seen that 
the contribution of the a term amounts to corrections 
of second-order O ((a — l)K r ), which we shall disregard. 
This is equivalent to set a = 1 in the above expressions. 
Hence, one obtains 



log £ 



/ 



K, 



Co 



d6 



' oa 7 



1 gc, 

Co 



« + "a" 



(47) 



Solving for £, yields 



c = Coe^+* 



(48) 



Obviously, setting K r = gives £ = £o, that is, one 
recovers the Schzarschild isotropic coordinates. One can 
now read the (isotropic) spatial component of the metric 
r~ 2 gee = r^sin^^g^ = g rr : 



5« = A(0 = e ^ 1 + 



M 
Wo 



(49) 



Recall that the Lorentz transformation maps r — > £, 
while £ is just a convenient integration variable. Hence, 
one must now write £o in terms of £. For this, one should 
invert the relation 



Since the perturbation is assumed to be small, one can 
consider that both K, K r <C r a . Hence, one expects 
£ also to be a small perturbation to the unperturbed 
isotropic coordinate £o! indeed, this is obtained setting 
K r = 0, from which follows 



K r 



(50) 



2?o 



and hence, to first order in K r 



and 



dr 



Vr 2 



2Mr 



Co 



(44) 



(51) 



Inserting this expression into the above equation and 
expanding to first order in U = M/£, one gets 



Co 1 



M 



Co 



2Co, 



2M \ 

- M 

r ) 



(45) 



Hence, to solve the perturbed K r ^ case, one per- 
forms the coordinate transformation r — > £o, obtaining 



= e «o «)+* i + 



M 



2Co(0 



~ 1 + 2 1 



(52) 



Through a similar procedure, the g tt component can 
be found: 



6 



1 = 



= _ e 2^ - 1 (l+^)" 



(53) 



3A- + Jf r \ rr2 



M 



where, as before, we have taken a = 1 (so that -ftT, 
have dimension of length). Both g tt and reduce to the 
usual isotropic metric components in the limit K, K r — > 
0. 

To read the PPN parameters one must now transform 
to a quasi-cartesian referential, where the metric reads 



ntt = -l + 2U-2f3U 2 
THj = (l + 2 7 [/)<% . 



(54) 



This is achieved by a (spatial) coordinate change £ — > £' 
so that £ = (1 — K/M)£'. This transforms the metric to 



f/tt = fltt 



-1 + 2 2- 



M 



M 



c -i + 2C/-|i-^i^:)[/ 2 



M 



9£ 
0f 



fftt 



1 + 2 1- 



Hence, one obtains 

~ 1 
7 ~ 1 

Inverting, one finds that 
K 



2K r 



M 



M ) (l-K/M)£>_ 
(55) 



K + K r 

M 
K + 2K r 
M 



(56) 



K r 

— ~ /? — 7 . 



(57) 



A drawback of these results is the dependence of these 
PPN parameters on K and K r , which are integration con- 
stant (more precisely, K r is defined as proportional to K, 
which is free valued). Therefore, these do not depend on 
the physical parameters associated with the breaking of 
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FIG. 1: Allowed values for K and K r , derived from \f3 — 1| < 
6 x 1CT 4 , 7 = 1 + (2.1 ± 2.3) x 10~ 5 and the assumption 
\1-K r /K\ <0.1. 



Lorentz invariance. This reflects the linearization of the 
Einstein's equation that was used in order to obtain the 
radially symmetric Birkhoff metric. In any case, one can 
conclude that the effect of temporal/radial LSB mani- 
fests itself linearly on the PPN parameters (3 and 7. 

The current bounds, arising from the Nordvedt effect, 
|/3 — 1| < 6 x 10~ 4 [18] and the Cassini-Huygens ex- 
periment, 7 = 1 + (2.1 ± 2.3) x 1(T 5 [19], can be used 
to constraint the parameter space (K,K r ). Taking the 
Sun's geometrical mass M = 1.5 km, one obtains the 
constraints 



\K- 
K A 



K r 



2K r 

Also, by definition 

K r = 



< 0.9 m , 

= (-3.1 + 3.4) x 1CT 2 



in 



lA-a 



(A + B) 



K 



(58) 



(59) 



with a ~ 1, A, B <C 1. Hence, one expects deviations of 
K r from K to be small. Thus, considering for instance 
the constraint |1 — K r /K\ < 0.1, the resulting range of 
allowed values for these parameters is depicted in Fig- 
ures 1 and 2. Notice that considering K <~ K r immedi- 
ately yields K ~ (—1 ± 1.1) x 10~ 2 , indicating that the 
perturbation has a very short range. Indeed, since this 
lies well inside the Sun, one should work with the inte- 
rior Scharzschild solution, and this value merely indicates 
that the effect outside the Sun is negligible. 

Returning to Eqs. (41) one notices that, in the limit 
M -> 0, one gets g tt = ~e 2K ^ and g rr = ~e 2Kr ^. This 
allows establishing an analogy with Rosen's bimetric the- 
ory, by interpreting these changes as due to a background 
metric r/^ u [14]. Notice that, in the absence of a cen- 
tral mass, the vector field no longer rolls to a radial vev, 
since this spatial symmetry is inherited from the van- 
ishing covariant derivative, that is, from the presence of 



7 



K r (m) P>»£ temporal and a spatial coordinate, assumed to lie on the 
K (m) x-axis, that is, = K _1 (a, 6, 0, 0). One calculates the 
metric perturbations to the flat Minkowsky metric. 



-0.005 



To first order in h^, one has 



-0.015 



(63) 



Roi 



(V 2 /i i - hko,ik) 



/ -0.025 L 

FIG. 2: Detail of Fig. 1, showing only the allowed region. 

the central mass. With this in mind, one can pursue the 
analogy with Rosen's bimetric theory and read the PPN 
parameter a 2 [15]: 

f]nv = diag (-CQ \ C7 1 , cr 1 , of 1 ) = ( M — o) ; ( 60 ) 
so that 



Rij — 2 h()0,ij ~t~ hkk,ij h>ki,kj h>kj,ki*) j 

where time derivatives were neglected, since one assumes 
that dCc. 

The axial LSB breaks the radial symmetry, and one 
must consider the Einstein equations: G M „ = T^+Tb^u, 
where Tb^v is the stress-energy tensor for the bumblebee 
field, 



Tb^v 



l -b a bf j R a[j9tlv - b^b a R m - b v b a R a 



(64) 



Its trace clearly vanishes and hence, from the trace of the 
Einstein equations one obtains 



a 2 = 



£o 

Cl 

2{K r - K) 



x = e 2(K r -K)/i _ 1 



a 



(A + B)- . 



(61) 



which exhibits a radial dependency. Assuming a ~ 1 and 
considering the spin precession constraint arising from 
solar to ecliptic alignment measurements [18], one has 
|a 2 | = \A + B\/r < 4 x 1CT 7 , which implies \A + B\ < 
4 x KT 7 r = 2.78 x 10 2 m, where r = 6.96 x 10 s m is 
the radius of the Sun. This also reflects the smallness of 
the expected effect on the metric. 

Since there is no explicit Lorentz breaking, the speed 
of light remains equal to c. The speed of gravitational 
waves, on the other hand, is shifted by an amount 



'£1 = e 2(K-K r )/(, 
Co 



A + B 



(62) 



and hence acquires a radial dependence. 

Let us close remarking that a direct comparison of the 
results of this section with the ones of the purely radial 
case is unfeasible. Indeed, the radial LSB effects are ex- 
act, while the radial/ temporal results are not. Therefore, 
one cannot directly consider the ^4 — > limit for compar- 
ison of the results obtained. 



IV. AXIAL/TEMPORAL LSB 

One now treat the anisotropic LSB case. For defini- 
tiveness, one assumes that the bumblebee field is stabi- 
lized at its vacuum expectation value and exhibits both a 



I)j„ + —g^T + Tb^v 



(65) 



First one calculates /loo to first order of the potential 
U\ for that one writes 



R 00 = -- (a 2 R m + b 2 R n + 2abR w )-2a (aR 00 + bR w ) 



Since 



and hence 



5o2\_o, b 2 
2 I ' 2 



1 + ^1 V 2 h 00 - ^00,11 = 



One can rewrite this as 



2 + 5a 2 - b 2 



2 + 5a 2 

This equation admits the solution 



(66) 



-Roo — —^ 2 hoo j Rm — , Rij — -hoo,ij , (67) 
one has 



i?oo = -2(« 2i? oo + & 2 i?ii)-2a 2 .Roo , (68) 



(69) 



/lOO.ll + ^00,22 + = . (70) 
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h 00 {x,y,z) = 



2M 



\jc\x 2 + y 2 + z 2 



(71) 



where c\ = (2 + 5a 2 )/(2 + 5a 2 - b 2 ). 

One now computes the components ha (i ^ 1) using 
this result: 



a, 



-2(2 + 5a 2 -6 2 )a 2 + 2(l + 2a 2 )6 2 



(2 + 5a 2 )(-2 + 6 2 ) 



, 1 + cd 



(78) 

for free ctj . In order to match with the unperturbed case 
r 2 = XiX 1 one chooses aj = 1 to get 



R H = - (a 2 R oa + b 2 Rn + 2abR Q1 ) , (72) 



hence 



ai = 1 



(2 - a 2 )b 2 



(2 + 5a 2 ) (2 -b 2 ) 



oti = otj = \ . (79) 



Hence, 



#oo 
Roi 

Rii 



00 



(V 2 /lii — hoo,ii + hu^ij 



hu(x,y,z) = h o(aix,y,z) 
where it has been defined 



2M 



yj c\x 2 + y 2 + z 2 



, (80) 



Ru — 2^°°. 11 ' 
and therefore (for j ^ i, j ^ 1) 



a 

y 



V 2 h 



0(1 



"2 I 2 2 /l "' 11 



One rewrites this equation as 



(73) 



(74) 



2 _ 2 2 _ 2(2 + 5a 2 - 26 2 )- 4a 2 6 2 
° 2 ~ aiC ° ~ (2 + 5a 2 - b 2 ){2 - b 2 ) 

We now compute the hu component: 



leading to 



(81) 



#n - 2 ( a2Roa + b2Rll "> ~ 2l>2Rl1 ' (82) 



(2 - b 2 ) h iiA1 + 4h iitii + 2h 



(a 2 - b 2 ) /ioo, n + (2 + a 2 ) h o,u + a 2 h Q0 



(75) 



33 



As an Ansatz for the solution, one takes hu(x,y,z) = 
hoo(aix,a2y,a3z). Notice that 



mjj = a 2 h Q0 j j (a 1 x,a 2 y,a 3 z) ~ Ojhoojj (x, y, z) (76) 



h 

and that 



and 



Roo 



1 + — j #ii — ~2~# 00 ' 



-V 2 h 



oo , 



(83) 



(84) 



(V 2 /m — ft.oo,n + hkk,u — 2/111,11) ■ 



As before, one can combine these equations into 



^00. 
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-2Mc, 



-2clx 2 +y 2 + z 2 
° (c 2 x 2 + y 2 + z 2 f 



l 00, 33 



-2x 2 + clx 2 + xl 
-2M 1 1 

(c 2 x 2 + y 2 + z 2 f 



(77) 



Substituting in Eq. (76) one obtains, after some calcula- 
tion (see Appendix I), the coefficients ac 



2 + a 2 + 36 2 

"11,22 + fill, 33 = 2 _[_ 352 "00,11 

a 2 

+ 2 + 352 (^00,22 + ^00,33) - 2/l22,ll • (85) 

Introducing the expression for /i o and /i 2 2 in the r.h.s. 
term, one gets 



hu,22 + /lll,33 = 2M 



« 2 (c^l) , 2 



2c\x 2 -y 2 - z 2 



2c: 



2 2c 2 x 2 -y 2 -z 2 



2 + 36 2 V (c 2 a; 2 + y 2 + z 2 ) 5 /2 2 (^2 + y2 + z 2)5/ 



(86) 
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It is clear that the solution is a linear combination of hoo 
and h 22 . Indeed, one has 



(87) 

+ c 2 , ) hoo(x,y,z) + 2c\h 22 {x, y,z) . 



hu(x,y,z) 

a 2 (c 2 o 1) 
2 + 36 2 



One now searches for the off-diagonal component h w . 
The component 



R w = -(a 2 + b 2 )R w -ab(R ao + R 11 ) , 
leads to 



(l + a 2 + b 2 )R w = -ab(R Q0 + Rn) , 



and 



Ron — — -^V hoo 



(88) 



(89) 



(90) 



-Rio 



l 



l 



(V 2 ft.i - ftio.n) , 

(V 2 /in - /ioo, ii + hkk,n — 2/iii,ii) 



Repeating the above procedure, one obtains 



h 



ab 



1,22 T "01.33 



1 + a 2 + b 2 

X [ft-00,22 + /i00,33 + /ill, 22 + /ill, 33 — 2/l22,ll] 



(91) 



Writing /ioi = — a&/ (1 + a 2 + & 2 )(/ioo + /in + <5/ioi) 5 ° nc 
finds 



<5/i 



01,22 



<5/l, 
„2„.2 



01,33 



2/l22,ll 
2\ 



and hence 



Therefore, 



4Mcl(2c|x 2 - y 2 - z 2 ) 
(c 2 x 2 + y 2 + z 2 ) 5 / 2 



<Sftoi(a:,j/,z) = -2c\h 22 {x,y,z) 



(92) 



(93) 



/ioi = 



ab 



l + a 2 + b 2 



2 + 36 2 + 



*oo 



(94) 



Finally, one computes /too to second order. The results 
are developed in the Appendix II. It is found that one 
gets only a correction to the first order term h$ : 



too 



2c 2 [6 + 96 2 + (15 + 22b 2 )a 2 ]+ ( 
c 2 (6+ 15a 2 + 6 2 )(2 + 36 2 ) 



l 00 



(95) 



Since the LSB clearly turns the metric anisotropic, the 
usual PPN paramctrization cannot be straightforwardly 
used to ascertain its effects. This is so as the PPN 
formalism relies on a quasi-cartesian frame of reference 
which, by definition, has all diagonal metric components 
gu equal. One might argue that there is a transformation 
to such a isotropic frame, but the obtained PPN param- 
eters would undoubtedly be unphysical. However, one 
can extract some PPN-like parameters from the results. 
First, one notes that 



y/(l + 2/c)cc 2 + y 2 + z 2 
For hoo, one gets 



(1 - kcos 2 0) . (96) 



hoo 



~ ( c o - l)cos 2 6 
2M ^-2 '- , (97) 



where no r~ 2 correction appears. Thus, one cannot ob- 
tain the PPN parameter (3. However, since hu ^ h 22 = 
hss, the same reasoning allows us to compute two pa- 
rameters analogous to the 7 PPN parameter. Recalling 
that 



h n (x,y,z) = - 
+ 2c\h< 12 {x, y,z) , 



a 2 (c 2 p-l) . , 
2 + 3b 2 + 



hoo(x,y,z) 



(98) 



one gets (after neglecting the normalization with respect 
to hoo), 



71 = 1 + cos 2 8 

\a 2 (c 2 -l) 



x 



2 + 36 2 

b 2 2 
l + —cos 2 9 , 



(1 - c 2 ) + 2c\{\ - c 2 ) 



72 = 1 + (1 - c\)cos 2 Q ~ 1 



ab 



cos 2 6 



(99) 



As expected, the effect of the x-axis LSB is mostly felt 
on the /in component. Direct comparison with the PPN 
parameter 7 is troublesome, given that the present case 
is obviously anisotropic. Hence, no clear connection can 
be derived to link 7 with 71 and 72. 

However, the measured value of 7 should be of the 
same order of magnitude as the average of 71 and 72, 
integrated over one orbit: 



10 



7 -l^( 7l+72 )-l^( COS 2 #)^(l- e 2 ) , (100) 

where e is the orbit eccentricity. For low values of e such 
as those found in the Solar System, one gets 7 ~ b 2 /8. 
Taking e ~ 0, the constraint 7 = 1 + (2.1 ± 2.3) x 10~ 5 
now yields |6| < 1.9 x 10~ 2 . 

As stated above, the standard PPN analysis fails in the 
present scenario, which is clearly anisotropic. A discus- 
sion involving anisotropy of inertia and its effect in the 
width of resonance lines can be found in Ref. [20] and 
references therein (see also Rcf. [21]). Presented as a test 
between Mach's principle and the equivalence principle, 
it relies on the hypothetical effect the proximity to the 
large mass of the galactic core could have on the proton's 
mass. By comparison, we note that a radial LSB with 
the galactic core as source would be perceived as an axial 
LSB in a region such as the Solar System. In Ref. [22] 
the bound Am/m < 3 x 10~ 22 can be found, m being 
the proton mass. Comparing with Eq. (97) gives 

£-(i-£)(.*-i)4s»xio-, (101) 

which yields |6| < 2.4 x 10 -11 , a much more stringent 
bound than the obtained above. 



V. CONCLUSIONS 

In this study we have obtained the solutions of a grav- 
ity model coupled to a vector field where Lorentz sym- 
metry is spontaneously broken. We have analyzed three 
cases: purely radial, temporal/radial and temporal/ axial 
LSB. 

In the first case, we have found a new black hole so- 
lution; we showed that, as in the standard scenario, this 
solution has a removable singularity at a horizon of ra- 
dius r s = {2Mtq L ) 1 ^ 1 ~ l , which is slightly perturbed 



with respect to the usual Scharzschild radius r s o = 2M. 
This horizon has an associated Hawking temperature of 
T = (2Mtq L )~ l Tq, and protects a real singularity at 
r = 0. Deviations from Kepler's law yield L < 2 x 10~ 9 . 

The temporal/radial case yields a slightly perturbed 
metric which leads to PPN parameters (3 « 1 — (K + 
K r )/M and 7 w 1 - (K + 2K r ) /M, directly proportional 
to the strength of the induced effect (given by K and 
K r ~ K). Unfortunately, since K and K r are integration 
constants, one cannot constrain the physical parameters 
from the observed limits on the PPN parameters. Also, 
analogously to Rosen's bimetric theory, one can obtain 
the PPN parameter 7 ~ (A+B)£, where £ is the distance 
to the central body and A and B rule the temporal and 
radial components of the vector field vev. 

In the temporal/axial case one gets, as expected, a 
breakdown of isotropy, and hence a standard PPN anal- 
ysis is not feasible. However, we have defined the 
direction-dependent "PPN" parameters 71 ~ b 2 cos 2 9/2 
and 72 ~ a 2 b 2 cos 2 9/A (considering the LSB occurs in 
the direction of xi). Naturally, 71 <C 72- A crude 
estimative of the measurable PPN parameter 7 yields 
7 ss b 2 (l — e 2 )/4, where e is the orbit's eccentricity. This 
final case requires further study, as its effects cannot be 
properly accounted for an isotropic formalism such as 
the parametrized post-Newtonian one. However, com- 
parison with experiments concerning the anisotropy of 
inertia yields \b\ < 2.4 x 10" 11 . 

Acknowledgments 

JP is sponsored by the Fundagao para a Ciencia 
e Tecnologia (Portuguese Agency) under the grant 
BD 6207/2001. 

VI. APPENDIX I 

Substituting Eqs. (76) in Eq. (76), one gets the fol- 
lowing relation: 



((2 - 6 2 )a 2 + (b 2 - a 2 )) ftoo.ii + (4a 2 - (2 + a 2 )) h 00 . u + (2a 2 - a 2 ) hoojj = , (102) 

and hence 

((2 - b 2 )a\ + (b 2 - a 2 )) c\ (-2c 2 x 2 + y 2 + z 2 ) + (4a 2 - (2 + a 2 )) (-2x\ + c 2 x 2 + x 2 ) 
+ (2a 2 - a 2 ) (-2x 2 + c 2 x 2 + x 2 M ) = . (103) 
For this equality to hold, the coefficients of each coordinate must vanish, leading to 

2(2 + 5a 2 )[-(2 + 5a 2 )(l + 2a\ - 2a 2 - a 2 ) + (-1 + 2a\ + a 2 (-4 + 5a 2 ) - 2a- - a 2 )b 2 ] = , 
-2(2 + 5a 2 )(2 + a 2 - 4a 2 + a 2 ) + [a 2 (-4 + 5a 2 ) + 2(1 + a\ - 4a 2 + a 2 )]b 2 = , 

-2(2 + 5a 2 )(-l + a\ + 2a 2 - 2a 2 ) + [a 2 (-4 + 5a 2 ) + 2(-2 + a 2 + 2a 2 - 2a 2 )]fe 2 = . (104) 
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These equations admit the solution 



-2(2 + 5a 2 - b 2 )a 2 + 2(1 + 2a 2 )& 2 2 1 + a 2 



ai = (2 + 5a 2 )(-2 + & 2 ) ' a< = -2- ' (1 ° 5) 



for a free ay. 



VII. APPENDIX II R 0i = R^) , since it contains no hoo term. The equation 



To obtain hoo to second order, one writes hoo — hoo + 
h$ and R 00 = R$+R$, R l3 = R^+R^f. Notice that 



for R$ is 



One obtains 



Hence 



and 



#oo =~\ (a 2 R$ + b 2 R ( ^ + 2abR w ) - 2a 2 R$ - 2abR w 



p(l) , p(2) 
"00 ' -"-oo 



„2 R (1) 

a R oa 



b 2 R[\ } + 2abR 10 ) - 2a 2 R 00 ) 



2abR 10 



+ \ (« 2 4o + b 2 R$ + 2abR w ) h$ + \ (a 2 R$ + b 2 R {2 ^ (-1 + h$) - 2a 2 R t 



2p(2) 
00 



(1 + 2a 2 ) + 2abR 10 



4o } + « 2 



5 , 
2 +h oo 



(i) 



i? (2) + - (- 
Hm + 2 I 



4 oo J 



-V J /i 



21,(2) 



-iv 2 ^ + 2ab (V 2 /iio - /iio,u) 



(i) 

00 



(l + 2a 2 ) 

i (V 2 hu - h$ >u + 2/122,11 - fcn,n) (-1 + 4o) 



2 



Dropping the (— 1 + h^) term yields 



It follows that 



2 



(106) 



(107) 



R$ = ((1 + 2a 2 )R^ + 2abR w ] h$ + ^R™ (-1 + h$) , (108) 



(109) 



(110) 



1 + 



bo 2 + 6 s 



' l 0041 ' V 1 ' 011 ) \ ""00,22 + "-do, 33 J 



~~2 (/in, 22 + /in, 33 + 2/122,11) 



(111) 



Taking the Ansatz = Ah^ , one finds 
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2MAc 2 (6 + 15a 2 + b 2 ) 2c 2 x 2 -y 



2 + 5a 2 



and 



(cgx 2 + y 2 + z 2 ) 



(5/2) 



-2M 



fe 2 [g 2 (c 2 - 1) + c 2 (2 + 3fr 2 )] 2c 2 x 2 - y 2 - z 2 

(2 + 5a 2 )(2 + 3fe 2 ) (c 2 :r 2 + y 2 + z 2 ) 5 / 2 ' 



(112) 



frV(c 2 -l) + c 2 (2 + 3fr 2 )] 
c 2 (6 + 15a 2 + 6 2 )(2 + 36 2 ) ' 

Thus, one obtains a small correction to Hqq , but no change in behavior. One can write 



(113) 



h -hW i ^ 2 > - 2c 2 ,[6 + 9b 2 + (15 + 22& 2 ) a 2 ] 



2l2 



a 2 6 : 



c 2 (6 + 15a 2 + 6 2 )(2 + 36 2 ) 



~ [ 1 - 



'00 



(114) 
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